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A brief report is presented on the effect of viscoelastic losses in a high density contrast sonic band- 
gap material of close-packed rubber spheres in air. The scattering properties of such a material are 
computed with an on-shell multiple scattering method, properties which are compared with the 
lossless case. The existence of an appreciable omnidirectional gap in the transmission spectrum, 
when losses are present, is also reported. 
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The problem of elastic wave propagation in inhomo- 
geneous media is of great importance in many branches 
of physics, mathematics and engineering. In particular, 
matters such as the localization of classical wavesEl and 
the formation of spectral gaps in periodic elastic compos- 
ites,a0 have drawn the attention of researchers over the 
last decade. In particular, phononic or sonic crystals are 
composite materials which consist of homogeneous par- 
ticles (solid or fluid inclusions the dimensions of which 
are large enough in order for a macroscopic description 
of their elastic properties to be valid) distributed pe- 
riodically in a host medium characterized by different 
mass density and Lame coefficients. With an appropriate 
choice of the parameters involved one may obtain sonic 
crystals with absolute frequency gaps (omnidirectional 
sonic gaps). 

Among the various methods available for the calcula- 
tion of the elastic properties of phononic crystals, the 
traditional band-structure methods mainly deal with pe- 
riodic, infinite, and nondissipative structures. However 
in an experiment, one deals with finite-size slabs and 
the measured quantities are, usually, the transmission 
and reflection coefficients. Apart from that, realistic 
structures are dispersive and have losses. This limita- 
tion of traditional band-structure calculations was no- 
ticed in a_theoretical study of colloidal crystals with ul- 
trasounds We remember that a usual band-structure 
calculation proceeds with a given wave vector, k, and 
computes the eigenfrequencies within a wide frequency 
range together with the corresponding eigenmodes. On 
the contrary, on-shell methods proceed differently: the 
frequency is fixed and one obtains the eigenmodes of 
the crystal for this frequency. These methods are ideal 
when one deals with dispersive materials (with or with- 
out losses). Moreover on-shell methods are computation- 
ally jpore efficient than traditional band-structure meth- 
ods, a Psarobas et al. have recently developed an on-shell 
method for phononic crystals ,Q which applies to systems 
which consist of nonoverlapping homogeneous spherical 
particles arranged periodically in a host medium charac- 
terized by different clastic coefficients. The method pro- 
vides the complex band structure of the infinite crystal 
associated with a given crystallographic plane; and also 
the transmission, reflection, and absorption coefficients 



of an elastic wave incident at any angle on a slab of the 
crystal, parallel to a given plane, of finite thickness. 

The present work introduces the effect of viscoelastic- 
ity in a sonic band-gap material. For this purpose, a 
binary system of close-packed rubber spheres in air is 
chosen. The viscoelastic response of the system is ac- 
counted for by means of the Kelvin- Voigt modelJ3 which 
is well-suited for materials and ultrasonic frequencies of 
major interest. The problem of acoustic- wave scattering 
by a single viscoelastic spherft-pf radius S has been ade- 
quately addressed in the pastil according to the Kelvin- 
Voigt viscoelastic model. In such a case the sphere is 
considered to be elastic with modified shear and com- 
pressional complex wavenumbers, the imaginary parts of 
which represent a measure of the loss. In particular, for 
an absorbing sphere in an inviscid fluid background, the 
complex compressional and shear wavenumbers are con- 
veniently defined as 
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where qi = w/cj refers to the fluid environment with u 
being the angular frequency and q the respective speed of 
sound. The real parts of the complex Lame parameters of 
the sphere, A s = \ se — i\ sv and p s = fi se —ip, sv , combined 
with the density p s yield the compressional and shear 
wave speeds respectively 
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The imaginary parts of the Lame parameters are con- 
nected to the viscous losses a + 2/3, and (3 of the sphere 
as follows: a — lj\ sv , f3 = lo[i sv . The problem aLelastic 
scattering by a solid sphere in an inviscid HuidEJ is de- 
scribed by the scattering transition matrix, the elements 
of which, in the angular momentum representation (I, m) 
(see Appendix), connect the spherical wave expansion co- 
efficientaj of the scattered field to those of the incident. 

Multiple scattering effects within planes of spheres, 
sonic crystals and slabs of the same are taken into account 
by the method described in Ref. @. This method views 
the crystal as a sequence of planes of spheres parallel to 
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a given surface: a crystallographic plane described by a 
two-dimensional (2D) lattice {R„}. The corresponding 
2D reciprocal lattice we denote by {g}. In the host re- 
gion between the n-th and the (n + l)-th planes, a Bloch 
wave solution for the displacement field (harmonic time 
dependence is assumed), corresponding to a given fre- 
quency to and a given reduced wave vector k|| within the 
surface Brillouin zone (SBZ) of the given surface, can be 
expanded into plane waves propagating (or decaying) to 
the left and to the right, as follows: 



u g „ exp 



iK- r (r-A„)}, (3) 



where 



K± = k||+g, ±[(c/ Q ) 2 -(k||+g) 2 ] 



2lV2 



(4) 



and A„ is a point between the nth and (n + l)th planes. 
We note that the separation of two successive planes of 
spheres need only be larger than the radii of the spheres 
of the two planes (see pg. 86 of Ref. ||). It should be 
also noted that, although both shear and compressional 
modes are considered within the spheres and enter the 
calculation through the scattering transition matrix (see 
Appendix) , only longitudinal waves exist in the host (air) 
region [see Eq. (0)]. As a consequence of that, in a binary 
composite of nonoverlapping solid spheres in a fluid host, 
such as the one being investigated in this report, there 
are no propagating shear waves if the solid component 
does not form a continuous network. 

A generalized, i.e. propagating or evanescent, Bloch 
wave satisfies the equation 
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where a 3 = A n+ i — A n and k = (kii,fc z (a;;k||)) is the 
Bloch wavevector. There are infinitely many such solu- 
tions for given k|| and u>, corresponding to different values 
of the z-component, fc z (w;k||), of the reduced wave vec- 
tor k, but in practice one needs to calculate only a finite 
number (a few tens at most) of these generalized Bloch 
waves. We have propagating waves [for these fc z (w;k|i) 
is real] which constitute the normal modes of the infinite 
crystal; and evanescent waves [for these k z (u)] k||) is imag- 
inary] which do not represent real waves, but they are 
useful mathematical entities which enter into the eval- 
uation of the reflection and transmission coefficients of 
a wave, with the same lo and k|| , incident on a slab 
of the crystal parallel to the given surface. The trans- 
mission/reflection matrices for a slab which consists of 
a stack of layers of spheres with the same 2D periodic- 
ity parallel to a given plane of the crystal are obtained 
from the transmission/reflection matrices of the individ- 
ual layers in the manner described in Ref. ^. Knowing 
the transmission/reflection matrices for the slab we can 




FIG. 1: (a): Normalized total scattering cross section of a 
rubber sphere in air. The solid line corresponds to a sphere 
with no losses and the shaded curve to a sphere with losses. 
The two different viscous levels used here, for the sphere, are 
hardly distinguishable. The resonance with the arrow ex- 
tends twice as high, (b) : Normalized absorption cross section 
of a viscoelastic rubber sphere. The dotted (solid) line corre- 
sponds to the low (high) viscous level used here. 



readily obtain the transmission, reflection, and absorp- 
tion coefficients of a plane acoustic wave incident on the 
slab. 

The system which will be examined here is a crystal 
of rubber spheres in air. The physical parameters en- 
tering our calculations are taken from Ref. 0. In par- 
ticular, the mass density of air is p a i r = 1.2 kg/m 3 and 
its respective speed of sound c air — 334 m/s. The rub- 
ber spheres have a mass density p s = 1130 kg/m 3 and 
cis = 1400 m/s, Cts — 94 m/s are the compressional and 
shear speeds of sound, respectively. In addition, accord- 
ing to Ref. 0, three different viscoelastic levels are con- 
sidered for the rubber spheres, namely: lossless spheres 
(a = (3 — 0), a low viscous level (oti ow = 0.5 MPa/s, 
Plow — 0.01 MPa/s) and a high viscous level (cthigh — 
5 MPa/s, fthigh = 0.1 MPa/s). The viscoelastic proper- 
ties used in this study are typical values for commercial 
rubbers, the variety of which is quite extensive^ aad fre- 
quency dependent at high ultrasonic frequencies. £21 

Before dealing with the main object of this work, it 
will be useful to state briefly certain basic features 
the acoustic scattering problem by a viscoelastic sphere 
These features will be used extensively in what follows. 
These are presented in Fig. ^, obtained with an angular 
momentum cutoff, l m ax = 5, which yields an accuracy of 
better than 0.001% in the given frequency range. On first 
sight, one may observe the disappearance of the sharp 
modal resonances in the scattering cross section of the 
sphere (see Appendix), when viscoelasticity is turned on. 
Also, there is no significant difference between the two 
viscous levels (low and high) in scattering. On the con- 
trary absorption (see Appendix) seems to be different for 
the two cases. The low viscous level induces resonant ab- 
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FIG. 2: The sonic band structure at the center of the SBZ of 
a (111) surface of an fee crystal of close-packed lossless rubber 
spheres in air (a). The corresponding transmittance curve of 
a slab of 16 layers parallel to the same surface is given in (b). 
In (c) the same transmittance curve is presented but with 
spheres of the low viscous level, d is the distance between 
successive (111) planes of the fee crystal under consideration. 



sorption exactly at the resonant frequencies of the lossless 
case, while the higher viscous level washes everything out 
as if there was no inner resonant structure in the system. 

We next consider an fee crystal of close-packed (almost 
touching) rubber spheres in air. We view the crystal as 
a succession of planes of spheres parallel to the (111) fee 
surface. Fig. |(a) shows the frequency band structure 
normal to the fee (111) plane (k|| = 0) and the corre- 
sponding transmission spectrum for waves incident nor- 
mally on a slab of the crystal consisting of 16 layers of 
lossless spheres. The results are obtained with an angular 
momentum cutoff l m ax = 7 and 55 g vectors (the estab- 
lished convergence is within an accuracy of better than a 
tenth of a percent). One observes, besides a large Bragg 
gap extending from LuS/c a i r = 1.223 to LuS/c a i r = 2.065, 
a number of flat bands which derive from the interacting 
sharp resonant modes localized on the individual rubber 
spheres (see Fig. |l|). Because these bands are so narrow in 
the present case, they are hardly observable; except that 
they introduce small gaps, above and below the main gap, 
which result from the hybridization of these flat bands 
with the broad bands corresponding to nearly free propa- 
gating waves. These narrow gaps are seen more clearly in 
the transmission spectrum [Fig. g|(b)]. Within the main 
gap these flat bands manifest themselves as sharp peaks 
in the transmission spectrum. The long wavelength limit 
is represented by the linear segments of the dispersion 
curves of Fig. g(a), the slopes of which determine the 
propagation velocity of acoustic waves (q = 1.54 c a ir) 
in a corresponding effective medium. The oscillations in 
the transmission coefficient, over the allowed regions of 
frequency, are due to interference effects resulting from 
multiple reflection at the surfaces of the slab of the crystal 
(Fabry-Perot-type oscillations). When losses are present 
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FIG. 3: Absorbance and transmittance curves of slabs of the 
rubber sonic crystal described in Fig. ^|(a) consisting of 8 
[(a),(c)] and 32 [(b), (d)] planes of spheres, respectively. The 
black line (shaded curve) corresponds to the low (high) vis- 
cous level. 



in the system, there are no true propagating waves and 
the band structure of the infinite lossless crystal is not 
of any help, therefore the effect of the low viscous level 
is shown in the transmission spectrum [Fig. 0(c)]. As 
expected from the results of the single sphere, the sharp 
peaks and dips of the resonant states disappear and we 
obtain a "clean" sonic gap without any resonant modes 
within it. The existence of the frequency gap means that 
sound does not propagate through the crystal when its 
frequency lies within the gap (the intensity of the wave 
decays exponentially into the crystal for these frequen- 
cies), and if it cannot enter into the crystal, it cannot be 
absorbed either. This is shown in Fig. |3[ 

The close relation between absorption and transmis- 
sion is demonstrated in Fig. || for slabs consisting of 8 
and 32 (111) planes of spheres, and for normal incidence. 
The solid lines correspond to the lower viscous level and 
the shaded curves to the higher one. A relatively large 
transmission coefficient implies that a correspondingly 
large fraction of sound has gone through the slab, with a 
consequent higher probability of being absorbed. How- 
ever, since losses are due to the rubber spheres, absorp- 
tion mainly occurs about the frequency regions where 
the modes of the acoustic field are mostly localized in 
the spheres. This explains why, outside the gap regions, 
absorption takes place essentially about the frequencies 
of the flat bands. 

In Fig. Q the projection of the frequency band struc- 
ture on the SBZ of the (111) plane of the fee crystal along 
its symmetry lines is shown. This is obtained, for a given 
k|| , as follows: the regions of frequency over which there 
are no propagating states in the infinite crystal [the corre- 
sponding values of of all k z (u>, k|| ) are complex] are shown 
in white, against the shaded areas which correspond to 
regions over which propagating states do exist [for a given 
uj there is at least one solution corresponding to k z {u>, k||) 
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FIG. 4: Projection of the frequency band structure on the 
SBZ of the (111) surface of the fee sonic crystal described in 
the caption of Fig. 0. The white areas show the frequency gaps 
in the considered frequency region. The broken lines map the 
omnidirectional frequency gap in the transmission spectrum 
of a slab of the crystal, with losses, of finite thickness. Finally, 
the inset shows the SBZ of the fee (111) surface. 



real] . One clearly sees here how the resonances on spheres 
lead to narrow hybridization gaps above and below the 
main gap, and flat bands in the gap regions. When losses 
are present, the crystal under investigation exhibits an 
appreciable omnidirectional sonic transmission gap ex- 
tending from coS/cair = 1.595 up to 1.946. Finally, we 
note that the results obtained depend on the ratio Xq/S 
(Ao is the wavelength of the incident sound wave) and are 
therefore applicable over different ranges of frequencies, 
provided that the viscoelastic properties of the spheres 
do not vary significantly with the frequency. 



Acknowledgments 

This work has been supported by the Institute of Com- 
munication and Computer Systems (ICCS) of the Na- 
tional Technical University of Athens. Support from the 
University of Athens is also acknowledged. 



APPENDIX A 



Applying the proper boundary conditions at the iiu 
terface between the surrounding fluid and the sphere,E!J 
requiring the continuity of the radial component of the 
displacement field and the su rface traction, along with 



the requirement that there is no tangential component of 
the surface traction at the interface, we can determine 
the scattering matrix which connects the incident with 
the scattered field. Here for reasons of completeness, 
along the lines of the formalism established in Ref. [(], 
the nonzero elements of the T matrix for a solid sphere 
in a fluid host are 
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with Z\ = Sqi referring to the fluid and x v — Sq sll , with 
v = l,t to the sphere. The superscripts LL refer to the 
case of scattering in a fluid host, since incident and scat- 
tered waves are L-type compressional waves. The 3x3 
determinants Di and W^ L are given by 
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d 2 2 = zihj (zi), d 23 = 1(1 + 1) Jit (act), d 24: = xij[(xi), 
d 32 = 0, d 33 = [1(1 + 1) - 1 - x 2 t /2] ji(x t ) - xtjftxt), 
du = xij[(xi) -ji(xi), d 4 2 = -x%phl(zi)/(2p a ), 

(1*3 = 1(1 + 1) [Xtj'i(xt) - jl(xt)} , 

d 44 = [1(1 + 1) - x 2 J2] j t (xt) - 2xij[(xi), 

zd'l(zi), 4 = 0, 4 = -x 2 t pj l (z l )/(2 Ps ) . 

(A3) 
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j[ and hf denote the first derivatives of the spherical 
Bessel and Hankel functions, respectively. The T ma- 
trix, because of spherical symmetry is diagonal in I and 
independent of m. 

The exact form of the above T matrix allows us to 
compute the normalized total scattering cross section of 
an elastic sphere (scattering cross section over irS 2 ) in a 
fluid host, 
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In addition the normalized absorption cross section is 
given by 
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(A5) 
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